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Method of determining an illuminated surface 
TECHNICAL AREA 

5 The present invention relates to a method of determining a surface illuminated by 

incident light by recording the intensity of light reflected from the area in a first image 
thereof and recording the intensity of light reflected from the area in a second image 
thereof, complementary to the first image, taken with another angle of illumination. 

10 The invention is particularly but not exclusively applicable to paper surfaces intended for 
the application of print. 

THE PRIOR ART 

15 

SE 508 822 makes known a method and a device for measuring and quantifying surface 
defects, such as polishing roses that can occur in connection with the polishing of coated 
sheet metal items. In this method and device, at least two sub-images are recorded with 
at least one camera under illumination of the test surface with parallel light or light from 

20 a point source, whereby the angles of incidence of the light relative to the test surface 

and/or the camera are different during the recording of different sub- images, after which 
the sub-images are processed in at least one central unit. After this, one or several 
difference images of the sub-images are produced, and used to determine the degree of 
surface defects on the test surface. This known technique, however, provides no 

25 guidance in how the recorded sub-images can be used in order to determine the 
topography of the surface. 

DESCRIPTION OF THE INVENTION 

30 

An object of the present invention is to provide a photometric method of the type 
described in the introduction that can rapidly determine the topography of a surface. 
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According to an aspect of the invention, the intensity (that is, the power per unit area) 
only of diffusely reflected light is recorded in the two images, and a difference between 
the recorded intensities of the diffusely reflected light of the first and the second recorded 
5 images is determined, in order to obtain a representation of the gradient variations of the 
surface. 

If the difference is normalised by division by the sum of the intensities, a ratio is 
obtained that is essentially directly proportional to the local derivative of the surface. 



The derivative in turn is used to determine the height function of the surface. 

The insight that forms the basis of the invention is that the lightness of a topographic 
surface element depends both on its diffuse reflectance and on its angle relative to the 
15 illumination. If images of the surface are taken with different angles of illumination, 
these will differ due to the topography of the surface, but not due to differences in its 
diffuse reflectance. This can, according to the invention, be used in image processing 
operations that distinguish the topography from the reflectance. 

20 

BRIEF DESCRIPTION OF DRAWINGS 

The invention is described in more detail with reference to the attached drawings, in 
which FIG. 1 shows schematically an arrangement for recording images according to the 
25 invention; FIG. 2 shows a model corresponding to FIG. 1 that forms the basis for 
processing the recorded images; FIG. 3 shows in the form of a diagram a simplified 
example of processing a recorded image accenting to the invention; FIGS. 4 A and 4B 
show images of a deeply printed test surface recorded by illumination from the left and 
from the right of the arrangement according to FIG. 1; FIG. 5 shows the reflectance of 



30 the test surface according to FIG. 4; FIG. 6 shows the derivative of the test surface in 
FIG. 4 ; FIG. 7 shows the topography of the test surface in FIG. 4; FIG. 8 shows an 
image of the test surface according to FIG. 4 with contours representing -1 ^tm; FIGS. 
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9 A and 9B show at a higher scale a reflectance image and a topographic image, 
respectively, of a test surface furnished with printed points; and FIG. 10 shows profiles 
of a test surface measured mechanically and measured with an arrangement according to 
the invention. 



DESCRIPTION OF AN EMBODIMENT 

10 The principle of the invention is shown in the arrangement according to FIG. 1 . A test 
surface 1, which in the examples described is a paper surface with an area typically of 
5x5 mm, is illuminated by a first light source 2 and by a second light source 3 arranged 
at two mutually opposite directions. The light sources 2, 3 contain halogen lamps with 
illumination optics. A camera 4 of CCD type detects and records by a computer 5 the 

15 intensity of the reflected light. 

The computer 5 is preferably equipped with known hardware and software for image 
processing. The time required for analysis of an image with a resolution of 512x512 
pixels is currently approximately 10 seconds using a 400 MHz standard PC. The 
20 mathematical analysis has been carried out using the MATLAB® software. 

The invention is based on the detection of diffuse light. Specular reflections from the test 
surface can be eliminated in the example shown by means of mutually crossed polarizers 
6 and 7. In more detail, a polarizer 6 can be placed between the test surface 1 and each 
25 light source 2,3, while a polarizer 7 that is crossed with respect to the polarizer 6 can be 
placed between the test surface 1 and the camera 4, in such a way that the illuminating 
light is polarised parallel to the incident plan** and the reflected light is polarised at right 
angles to it. 



5 



30 



With reference to FIG. 2, the intensity of the incident light is proportional to cos(a), 
where a is the angle of incidence of the illuminating light to the surface 1. Lambert's law 
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is assumed to be valid for the diffusely spread light. According to this law, the radiance 
is equal in all directions. This means that the intensity detected by the camera is given by 

I = Io R cos(a) [1] 

where R is the reflectance and Io is the intensity measured when R = 1 and cos(a) = 1 . 
Scalar multiplication gives a value for cos(a) as 



10 cos(ar) = a • n/|n| = 



sin(^)— + cos(x) 

dx 



+ 1 





2 

+ 











[2] 



where a is the illumination vector [sin(y),0,-cos(y)] and n is the surface normal 

[df/dx^df/dy,-!] 



15 If two images, L and I2, are recorded with 72 = -yi, FIGS. 4A, 4B, the partial derivative 
df /dxcan be calculated from [1] and [2] as 



3f = 1 h~h 
ck tan y I x + I 2 



[3] 



20 This expression does not depend on the reflectance. An example of the derivative, 

calculated from the images in FIGS. 4A, 4B, is shown in FIG. 6, in which the derivative 
has been coded as a grey -scale image. 



In order to obtain the height function of the test surface, the derivative must be 
25 integrated. However, since the images contain noise, certain spatial frequencies must be 
integrated with caution. This is why the derivative should preferably be subjected to a 
Fourier transform and multiplied by what is known as a Wiener filter: 

H * = , ,2 "'' 1 W 

\H,\ +SNR(u,vy* 

30 which performs the integration with the suppression of spatial frequencies u and v, which 
have an expected low signal-to-noise ration, SNR. The frequencies Hi of the surface 
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include both the partial derivative (in the form of 2mu) and the light that is spread in the 
material. For more detailed description of a Wiener filter, refer to Pratt, W. K., (1978), 
Digital Image Processing, Wiley, New York, 378-387. The surface function, which is 
shown in FIG. 7, also coded as a grey-value image in which lower surface areas have a 
darker grey value than higher surface areas, is obtained as the inverse transform of the 
product. 



The local reflectance of the test surface, which provides information about the degree of 
covering of the print, is approximately obtained as the sum of the images, Ii and L, see 
10 FIG. 5. 



In order to facilitate understanding of the invention, a simplified one-dimensional 
"digital" observation of a typical image processing operation is shown in FIGS. 3 A-G. 

15 FIG. 3 A shows the test surface, the topography of which, f(x), is to be investigated. In 
this case the surface has a printed regular pattern. 

When the surface is illuminated with oblique illumination from the left, an intensity 
variation in the diffusely reflected light is obtained, according to FIG. 3B, as a result of 
20 variations both in the reflectance (the pattern) and in the topography. Compare also the 
equivalent image or graphical representation in the two-dimensional case according to 
FIG. 4A, in which variations in grey- value are equivalent to variations in intensity. 

When the surface is illuminated with oblique illumination from the right, a new intensity 
25 variation, h(x), in the diffusely reflected light is obtained in an equivalent manner, 

according to FIG. 3C. Compare also the equivalent image in the two-dimensional case 
according to FIG. 4B. 

If the difference between the intensities, Ii(x) - h(x), is calculated, a variation which 
30 accentuates the topographic variations is obtained, according to FIG. 3D (the variations 
in reflectance are partially, but not wholly, suppressed), that is, variations in the gradient 
of the surface. 
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If the sum of the intensities, I,(x) + h(x), is calculated, a variation that essentially 
depends only on variations in reflectance is obtained according to FIG. 3E, while the 
structural or topographical variations are suppressed. In other words, the distribution of 
5 colour on the surface is obtained, that is, the presence or absence of print. Compare also 
the equivalent image in the two-dimensional case according to FIG. 5. 

If the ratio (Ii(x) - b(x))/( I.(x) + h(x)) is calculated, that is, the normalised difference 
between the intensities, a variation is obtained according to FIG. 3F that essentially only 
10 depends on topographic variations, that is, variations in the gradient of the surface. 

The ratio is used to calculate the derivative of the surface according to FIG. 3F as 

r( , __L_ /,(*Wi(*) 
AW ~tany 7,(x) + / 2 (x) 

where y = the angle of incidence of the illumination, as previously. Compare also the 
15 equivalent image in the two-dimensional case according to FIG. 6 In the two- 
dimensional case the derivative will in an equivalent way become 

1 /,(x,.y)-/ 2 (s,J') 
AU,j;) *tany I x (x,y) + I 2 {x,y) 

If the derivative is integrated, preferably with the simultaneous suppression of noise as 
described above, the topography is obtained according to FIG. 3G. Compare also the 
20 equivalent image in the two-dimensional case according to FIG. 7 

As the previous description has made clear, in addition to the purely topographic 
determination (FIG. 7) of a surface, the invention can also be used for the simultaneous 
determination of the reflectance of the surfac. (FIG. 5) in the same co-ordinates. In this 
25 way, interesting relationships between surface structure and the transfer of print can be 
studied in detail. In FIG. 8, contours equivalent to a depth of -1 urn from a sliding 
reference level have been added to the reflectance image from FIG. 5 by a thresholding 
operation in the image processing computer 5, which explains why printed points are 
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missing from regions of the printed area. In a similar manner, an examination has been 
made on the test surface according to FIGS. 9 A and 9B whether a particular depth of the 
depressions in the surface (for example, dark regions in the upper left-hand corner of the 
topographical map FIG. 9B) can correspond to failed print transfer (missing print points 
5 in FIG. 9A in those regions that are darkest in FIG. 9B). This can be used in printing 
technology as a prediction about in which regions missing printed points can be expected. 



In this respect it has become evident that so-called straight thesholding of the topography 
does not work so well. On the other hand, if a high-pass filter is applied to the 

10 topographical map such that long wavelength information is suppressed, and then apply a 
threshold level of -1 jim, that is to say, in practice thresholding relative to a sliding 
reference level, then the areas that have a high probability for missing print transfer are 
marked, see FIG. 8. It is possible to learn from this more about how surface rawness 
should be measured in a manner that is relevant for printability. The method has also 

15 given interesting results for full-tone areas printed flexographic printing (not shown). 

It is not necessary that the two images be recorded at different times. For example, the 
first image can be recorded with the arrangement described in FIG. 1 in a first 
wavelength region and the second image can be simultaneously recorded from the same 
20 camera point in a second wavelength region, complementary or distinct from the first 
wavelength region (not shown), if the two illuminations use distinct wavelength regions. 
In this way the possibility or recording processes on the test surface, for example, a 
region of a paper pathway during production, which is in motion. 

25 Analyses according to the invention of test pieces of LWC paper have shown a high 

correlation, r 2 =0.95, between profiles determined according to the invention and profiles 
determined according to conventional optical and mechanical methods of measuring 
profiles. In the diagram shown in FIG. 10, the full curve shows the profile determined 
according to the invention, while the dashed curve shows the same profile of the same 

30 paper strip determined by a mechanical contact method of measuring. 



